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Abstract 

This  paper  studies  the  problem  of  transistor  sizing  of 
CMOS  circuits  optimized  for  energy-delay  efficiency,  i.e., 
for  optimal  Etn  where  E  is  the  energy  consumption  and 
t  is  the  delay  of  the  circuit,  while  n  is  a  fixed  positive  op¬ 
timization  index  that  reflects  the  chosen  trade-off  between 
energy  and  delay. 

We  propose  a  set  of  analytical  formulas  that  closely  ap¬ 
proximate  the  optimal  transistor  sizes.  We  then  study  an 
efficient  iteration  procedure  that  can  further  improve  the 
original  analytical  solution.  Based  on  these  results,  we 
introduce  a  novel  transistor  sizing  algorithm  for  energy- 
delay  efficiency. 

1.  Introduction 

The  rapidly  increasing  complexity  of  VLSI  systems 
has  made  it  necessary  to  pay  ever  more  attention  to  de¬ 
sign  issues  that  affect  energy  consumption.  One  of  the 
original  motivations  for  CMOS  technology  was  its  low  en¬ 
ergy  consumption,  and  today,  there  are  still  no  alternatives 
that  approach  it  in  energy  efficiency.  Nevertheless,  energy 
consumption  is  more  and  more  often  the  factor  that  limits 
the  performance  of  contemporary  CMOS  systems. 

In  order  to  compare  designs  that  run  at  different  speeds 
and  consume  different  amounts  of  energy,  we  have  to 
combine  the  energy,  E,  and  the  delay,  t,  into  a  single  met¬ 
ric  S.  The  authors  have  previously  proposed  <5  =  Et 2 
as  an  energy-delay-efficiency  metric  for  VLSI  computa¬ 
tion  [1,  2,  17],  The  main  reason  for  choosing  this  met¬ 
ric  over  others  is  that  S  is  to  first  order  constant  when  we 
vary  the  supply  voltage  of  a  CMOS  system:  the  delay  falls 
roughly  linearly  with  supply  voltage,  and  the  energy  con¬ 
sumption  increases  roughly  quadratically;  therefore,  Et2 
stays  roughly  constant.  Hence,  the  S  metric  allows  the  de¬ 
signer  to  factor  “runtime”  voltage  scaling  out  of  consider¬ 
ation.  The  authors  have  argued  that,  owing  to  its  voltage 
independence,  the  5  metric  is  superior  to  other  efficiency 
metrics  found  in  the  literature,  such  as  E  or  Et  [3], 

In  practice,  we  can  achieve  any  desired  target  speed  or 
target  energy  consumption  by  adjusting  the  supply  volt¬ 
age.  If  we  desire  to  change  to  a  particular  delay  target  t, 
we  adjust  the  voltage  to  meet  it,  and  a  circuit  optimized 


for  6  would  have  the  best  E  for  that  t.  Likewise,  we  may 
choose  an  energy  target  E  and  get  the  best  t  instead. 

The  Et 2  metric  is  a  special  case  of  a  wider  class  of 
energy  metrics,  which  includes  E  and  Et,  among  others. 
The  authors  have  shown  that  a  metric  of  the  more  general 
form  Etn  for  n  >  0  characterizes  any  feasible  trade-off, 
not  only  the  trade-off  through  voltage  scaling,  between  the 
energy  and  the  delay  of  a  computation  [4],  For  example, 
any  problem  of  minimizing  the  energy  of  a  circuit  for  a 
given  target  delay  can  be  restated  as  minimizing  Etn  for 
a  certain  n.  We  call  n  the  energy-delay  efficiency  index. 

In  this  paper,  we  study  the  problem  of  transistor  siz¬ 
ing  for  energy-delay  efficient  circuits.  Given  a  transistor 
netlist  where  each  transistor  i  has  width  Wi  and  length  A 
transistor  sizing  finds  the  values  of  Wi  and  1-,  that  optimize 
the  target  function — in  our  case  Etn.  While  it  is  true  that 
most  layout  systems  demand  that  transistor  sizes  be  quan¬ 
tized  to  some  grid,  we  ignore  this  constraint. 

Also,  we  can  remove  the  /;s  from  consideration  since 
there  is  usually  no  reason  to  set  the  lengths  of  transis¬ 
tors  in  a  digital  circuit  to  anything  other  than  the  mini¬ 
mum  allowed  by  the  fabrication  technology:  increasing 
the  length  increases  both  the  resistance  and  the  capaci¬ 
tance  and  hence  worsens  both  the  energy  and  delay. 

The  sized  transistors  of  a  circuit  are  connected  to  each 
other  through  wires.  The  capacitance  of  these  wires  leads 
to  additional  energy  and  delay.  (We  ignore  wire  resis¬ 
tance  in  this  paper.)  For  delay-only  optimization,  which 
can  be  phrased  as  the  minimization  of  the  metric  Etn  for 
very  large  n,  the  wire  capacitance  can  be  overcome  by  in¬ 
creasing  transistor  sizes  where  appropriate.  Conversely, 
for  energy-only  optimization,  when  n  =  0,  the  transistor 
widths  can  be  chosen  to  be  minimum  size,  independently 
of  the  wire  capacitance.  In  contrast  to  these  special  cases, 
for  n  small  but  nonzero,  wire  capacitance  cannot  be  ig¬ 
nored  or  overcome  in  a  straightforward  way,  and  the  opti¬ 
mal  transistor  sizes  depend  strongly  on  this  capacitance. 

In  this  paper,  we  propose  an  analytic  formula  for  tran¬ 
sistor  sizing.  If  the  approximate  solution  is  acceptable  for 
the  given  application,  the  formula  can  be  used  as  is  (no 
numerical  optimization  is  then  needed);  however,  if  more 
accuracy  is  required,  the  formula  can  be  used  to  provide 
a  good  starting  point  for  numerical  optimization.  Later  in 
the  paper,  we  propose  an  efficient  iteration  procedure  that 
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can  further  improve  the  accuracy  of  the  original  analyti¬ 
cal  solution.  Based  on  these  results,  we  introduce  a  novel 
transistor-sizing  algorithm  for  energy-delay  efficient  cir¬ 
cuits. 

The  proofs  of  properties  and  theorems  have  been  omit¬ 
ted  owing  to  space  limitations.  They  can  be  found  in  the 
first  author’s  Ph.D.  dissertation  [16]. 

2.  Previous  Work 

Classical  numerical  methods,  such  as  the  conju¬ 
gate  gradient  descent  method,  have  been  applied  to  the 
transistor-sizing  problem:  there  exist  several  transistor 
sizing  programs  that  minimize  power  consumption  while 
maintaining  performance  specifications  [5,  6,  7].  More 
recently,  several  specialized  numerical  techniques  have 
been  proposed  [8,  9,  10],  On  the  analytical  side,  Cong 
and  Koh  have  studied  the  related  problem  of  simulta¬ 
neous  gate  and  wire  optimization  for  optimal  delay  and 
power  [13],  Cong  and  Koh’s  solution  space  and  opti¬ 
mization  metric  are  different  from  what  we  shall  see  in 
the  present  paper.  A  different  analytic  approach  to  the 
transistor  sizing  problem,  for  the  performance  metric  Et , 
is  given  by  Hu  [11]  and  another  by  Horowitz,  Indermaur, 
and  Gonzalez  [12].  Both  Hu  and  Horowitz  et  al.  present 
qualitative  results;  they  only  analyze  basic  inverter  gates. 
To  the  best  of  the  authors’  knowledge,  the  present  paper  is 
the  first  one  that  goes  beyond  such  a  qualitative  approach, 
both  in  terms  of  the  generality  of  the  optimization  metric 
and  in  terms  of  the  generality  of  the  considered  circuits. 

3.  Ei”-optimal  circuits 

Let  t  be  the  cycle  time  of  the  critical  cycle  of  the  circuit 
whose  transistor  sizes  we  wish  to  optimize.  We  assume 
that  the  circuit  is  designed  so  that  all  cycles  are  critical; 
this  is  true  in  many  well  designed  circuits,  and  it  is  true 
for  any  optimally  sized  circuit  in  the  absence  of  additional 
constraints  on  transistor  sizes  (such  as  minimum-size  con¬ 
straints  or  slew-rate  constraints).  Let  E  be  the  energy  con¬ 
sumption  of  the  critical  cycle.  Let  us  further  assume  that 
E  is  a  constant  proportion  of  the  total  energy  consump¬ 
tion;  in  this  case,  optimizing  the  energy  E  of  the  critical 
cycle  optimizes  the  total  energy  of  the  circuit,  and  vice- 
versa. 

Using  the  r-model  [14,  15],  we  can  write  the  energy  as 

m— 1 

E  =  ^(Wni+Wpi+Pi),  (1) 

i= 0 


and  the  delay  as 


^  knifi+l  (^n(«+l)  1)  Pi+1 ) 

i= 0 
m— 1 


Wni 


_l_  1  (tfn(i+l)  +  wp(i+l)  +  Pi+ 1)  ^ 


i= 0 


W. 


pi 


where  wni  and  wpi  are  the  nFET  and  pFET  (nMOS  and 
pMOS  transistor)  widths  of  logic  gate  i;  kpj  >  0  are 
the  numbers  of  nFETs  and  pFETs  in  series  in  logic  gate 
i;  p.j+ 1  >  0  represents  the  wire  parasitics  at  the  output  of 
logic  gate  i;  /j+i  >  0  is  the  fanout  of  logic  gate  i;  p  is  the 
ratio  of  electron  mobility  to  hole  mobility;  m  is  the  length 
of  the  cycle,  and  i  £  0 ..m  —  1  with  all  indices  modulo  m. 
In  writing  Equations  1  and  2  we  have  made  several  simpli¬ 
fying  assumptions.  We  ignored  the  energy  consumption 
due  to  short-circuit  and  leakage  currents.  Furthermore,  we 
have  constrained  all  devices  in  a  series  transistor  network 
to  have  the  same  width.  Finally,  we  have  ignored  the  wire 
RC  and  time-of-flight  delays. 


4.  Properties  of  transistor  sizes  in 
.Et” -optimal  circuits 


Property  1  If  Wi  are  the  values  that  minimize  Etn  for  a 
given  set  of  wire  parasitics  Pi  and  gate  topologies  hi,  then 
awi,  a  >  0  are  the  values  that  minimize  Etn  for  the  set 
of  wire  parasitics  api  and  gate  topologies  hi- 
Property  2  If  Wi  are  the  values  that  minimize  Etn  for  a 
given  set  of  wire  parasitics  Pi  and  gate  topologies  fcj,  then 
Wi  also  minimize  Etn  for  the  set  of  wire  parasitics  pi  and 
gate  topologies  aki. 

If  we  ignore  special  constraints  on  transistor  sizes,  such 
as  minimum-size  and  minimum-slew-rate  constraints, 
and  if  we  further  assume  that  every  transition  on  every 
circuit  node  matters  to  the  circuit’s  overall  speed  (this 
last  assumption  is  especially  relevant  in  asynchronous  cir¬ 
cuits),  then  we  can  show  that,  when  a  system  is  optimized 
for  Etn,  the  widths  of  the  nFETs  and  pFETs  of  each  gate 
i  are  related  as  follows  [16]: 


tUpi  —  tVpi 


(3) 


Equation  3  is  a  local  relationship;  it  does  not  depend  on  ei¬ 
ther  E,  t  or  n.  Equation  3  allows  us  to  eliminate  either  the 
nFETs  or  the  pFETs  from  the  transistor-sizing  problem. 
In  particular,  with  the  notation 


Wi  —  Wni  Wpi  —  Wni  I  1  "I" 


pkp 


and 


ki  —  fi+lkni  I  1  + 


pkpi 


(4) 

(5) 


by  eliminating  the  pFET  sizes  from  Equations  1  and  2,  we 
get 

m  —  1 

E='^2(wi+pi)  (6) 

i= 0 


and 

t=y  kiWi+1+Pi+1.  (?) 

U  Wi 

We  shall  use  these  simpler  formulas  in  the  expressions  for 
Etn. 


5.  Preliminaries  for  £’tn-optimal  transistor 
sizing 

We  formalize  the  sizing  problem  of  a  transistor  netlist 
for  minimal  Etn  as  the  minimization,  over  the  WiS,  of  Etn 
where  E  and  t  are  given  by  Equations  6  and  7. 

(TO  — 1  \  /TO  — 1  \  11 


6.  Etn -optimal  transistor  sizes 

So  far  we  have  explored  some  general  properties  of 
transistor  sizes  for  circuits  optimized  for  Etn.  Based  on 
these  properties,  we  now  develop  a  simple  analytical  for¬ 
mula  that  approximates  the  transistor  sizes  of  an  Etn- 
optimal  circuit. 

We  start  by  finding  an  approximate  formula  for  the 
transistor  sizes  Wi  that  optimize  Etn,  in  Equation  8,  when 
Vi  :  ki  =  k.  We  then  extend  this  formula  to  the  case  when 
the  kiS  are  no  longer  equal  to  each  other. 


Note  that  Equation  8  holds  not  only  for  a  ring,  but  also 
for  a  chain  of  gates,  as  long  as  the  widths  and  parasitics  for 
the  input  of  the  chain  are  equal  to  the  widths  and  parasitics 
for  the  output  of  the  chain  (since  in  this  case  the  E  and  t 
for  a  chain  have  the  same  form  as  the  ones  for  a  ring). 
This  is  an  important  observation,  as  it  makes  our  results 
for  transistor  sizing  applicable  both  to  latency  and  cycle¬ 
time  minimization. 

Et"  is  a  posynomial  function  of  the  transistor  widths. 
A  posynomial  in  variables  Wi  is  a  function  of  the  form 

oO  ol  Brn~1 

J2oci<qaiwo  wi  -wm-i  where  ct*  >  0.  A  posyno¬ 
mial  problem  is  the  minimization  of  one  posynomial  while 
simultaneously  satisfying  a  set  of  upper-bound  constraints 
on  other  posynomials.  With  the  substitution  Wi  =  eXi ,  any 
posynomial  can  be  transformed  into  a  convex  function; 
therefore  the  unique  optimum  of  Etn  is  achieved  when 

V*:^=0. 

This  implies  that  the  optimum  is  achieved  when 


Vi : 


ki- i  _  kj(wi+ 1  +pi+ 1) 
Wi- 1  wf 


V^TO-t  ,  Wj  +  Pi 

1  2-,i= o  u’i-i 
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nP ’ 
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where  P  =  E /t  is  the  power  consumption  of  the  chosen 
cycle.  If  Vi  :  Pi  =0  (no  wire  parasitics)  and  n  is  very 
large  (delay-only  optimization).  Equation  9  reduces  to 


k.w*+ 1  -u.  ,  Wi 

/i?  —  hi—  1  " 


Wj 


Wi- 1 


which  is  the  known  condition  of  equal  stage  delays  for 
delay-only  transistor  sizing  [14].  If  we  were  able  to  solve 
Equation  9  analytically  for  any  pL  s  and  fc,-s,  we  could 
compute  the  optimal  directly  and  our  transistor-sizing 
problem  would  be  solved.  Unfortunately,  this  is  not  the 
case.  We  can  compute  an  exact  analytical  solution  of 
Equation  9  only  for  a  restricted  class  of  p,s  and  ki  s  [16]. 
In  particular,  we  can  show  that  if  Vi  :  ki  =  k,  i.e.,  the  case 
of  homogeneous  circuits,  and  Vi  :  Pi  =  p  then 


Vi  :  Wi  =  np  Mp,  k  >  0.  (10) 


Equation  10  states  that  the  transistor  widths  Wi  of  a  homo¬ 
geneous  circuit  with  equal  wire  parasitics  p,  optimized  for 
Etn,  are  all  equal  to  np,  independently  of  k  [17,  18], 


6.1.  Homogeneous  Circuits 

For  the  case  when  Vi  :  ki  =  k,  we  propose  an  approxi¬ 
mate  solution  of  the  tfjS,  of  the  following  form: 


Wi  =  aiPi+i  +  a2PAvg  (1 1) 

where 

PAvg  =  —y^Pi  (12) 

m,  ^ 

and  cti  and  a2  are  constants  to  be  determined  later.  First, 
let  us  motivate  Equation  11.  Based  on  Property  2,  we 
know  that  finding  the  wj,s  when  Vi  :  ki  =  k  is  equiva¬ 
lent  to  finding  the  u^s  when  Vi  :  A:,-  =  1.  In  other  words, 
the  value  of  the  tfjS  is  independent  of  the  kp,  when  all  fcjS 
are  equal.  Conversely,  based  on  Property  1 ,  we  know  that 
the  wjjS  scale  linearly  with  the  pjS.  This  suggests  that  the 
should  not  have  terms  that  are  independent  of  the  p;s. 
Based  on  our  experience  of  sizing,  we  know  that — while 
the  transistor  sizes  of  gate  i  depend  mostly  on  pi+ 1 — the 
effect  of  a  particular  pi  gets  distributed  to  some  degree  to 
all  other  gates.  As  a  consequence,  we  would  like  Equa¬ 
tion  1 1  to  depend  linearly  on  both  pi+ 1  and  some  average 
of  all  other  pjS  and  one  such  choice  is  a.iPi+i  +  a2pAvg- 
We  use  the  arithmetic  mean  for  pAVg  since  the  p\  s  cor¬ 
respond  physically  to  wire  capacitances  that  are  manipu¬ 
lated  additively  both  in  terms  of  delay  and  energy.  With 
these  clarifications  in  mind,  we  state  the  following: 

Theorem  1  For  a  neighborhood  Vp  =  \p  —  r],p  +  rj\  of 
p  >  0,  rj  >  0,  the  values  ofa±  and  a2  that  minimize  Etn 
given  the  WiS  of  the  form  defined  by  Equation  11,  where 
Mi  :  Pi  £  Vp,  ki  =  k  >  0  and  r]  — >  0,  are 

t  l 

ai  =  i  |  2  to  and  «2  =  n  -  ,  2  m  ■ 

n  m— 1  n  m— 1 

If  the  problem  is  large,  i.e.,  m  — >  oo,  — ^  ta  1  =>• 

Q1  =  WTnj  and  «2  =  thus  sa  =  1  +  2n. 

What  is  particularly  surprising  about  Equation  1 1  is  that 
the  strength  of  a  given  gate  depends  far  more  strongly  (5  x 
for  Et 2  optimization)  on  the  average  parasitic  load  (a2  = 
5/3)  than  it  does  on  the  load  on  that  particular  gate  (cti  = 
1/3).  Furthermore,  lim„_j.o  cti  =  lim«-s.o  a2  =  0  =>•  Vi  : 
Wi  =  0  for  n  =  0  regardless  of  the  piS.  In  other  words,  for 
energy-only  optimization.  Equation  1 1  yields  minimum- 
size  transistors,  as  one  might  expect. 
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optimization  index  (n) 

Figure  1.  Accuracy  in  E,  t  and  Etn  of  Equation  11  with 
oti  and  «2  given  by  Theorem  1. 

Theorem  1  yields  the  optimal  values  of  ai  and  a2  in  a 
close  neighborhood  of  p,  or  equivalently  when  the  p,s  are 
close  to  each  other.  We  want  to  check  now  if  the  form  of 
the  WiS  given  by  Equation  1 1  and  Theorem  1  yields  a  prac¬ 
tical  approximation  of  the  ups  when  Vi  :  Ap  =  k  but  the 
PiS  are  no  longer  close  to  each  other.  We  use  a  numerical 
optimizer  to  compute  the  error  between  the  optimal  and 
the  predicted  Etn  for  a  given  n,  m  and  a  set  of  p.;s.  We 
varied  to  £  [2, 1000],  n  €  [1, 10]  and  used  three  differ¬ 
ent  distributions  (uniform,  uniform-squared,  and  uniform- 
cubed)  for  pi  (E  [1, 100].  The  observed  errors  are  practi¬ 
cally  independent  of  the  problem  size  to  and  the  distribu¬ 
tion  chosen  for  the  p, s;  the  errors  only  depend  on  n.  Fig¬ 
ure  1  shows  the  relative  error  in  E,  t  and  Etn  for  to  =  31, 
n  £  [1,10],  hi  =  1  and  pj  £  [1,100]  chosen  randomly 
through  a  uniform- squared  distribution.  The  average  error 
in  E  is  between  4.1%  and  5.5%,  the  average  error  in  t  is 
between  -3.0%  and  -0.3%,  and  the  average  error  in  Etn  is 
between  1.0%  and  1.7%. 

6.2.  Non-homogeneous  Circuits  (first  form) 

The  formula  resulting  from  Theorem  1  yields  excellent 
results  when  all  kj s  are  equal.  We  would  like  to  extend 
it  to  incorporate  the  case  when  the  Aps  are  no  longer  all 
equal.  To  do  this,  we  assume  that  the  cumulative  effect 
of  the  pjS  and  the  Aps  on  the  ups  can  be  viewed  as  the 
product  between  the  individual  effect  of  the  p,s  (wire  ca¬ 
pacitances)  on  the  ups  and  the  individual  effect  of  the  Aps 
(gate  topologies)  on  the  ups.  Hence,  we  propose  an  ap¬ 
proximate  solution  of  the  up  s  of  the  following  form: 

Wi  =  (aipi+i  +  aiPAvg )  ri(k0,  Ap, . . . ,  km- 1)  (13) 

where  op  and  n->  are  given  by  Theorem  1,  while  functions 
Vi  will  be  determined  later.  When  all  gates  are  identical, 
i.e.,  Vi  :  Ap  =  Ap  we  know  from  Equation  10  that  the  ups 
are  independent  of  the  Aps.  For  this  reason,  we  choose 
ri(ko,  Ap, . . . ,  km- 1)  such  that  Vfc  :  ri(k,  k,...,k)  =  1. 

Based  on  our  experience  on  delay-only  transistor  siz¬ 
ing,  we  know  that — while  the  transistor  sizes  of  gate  i 
depend  strongly  on  ki — the  effect  of  a  particular  ki  gets 
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Figure  2.  Accuracy  in  E,  t  and  Etn  of  Equation  13  with 
Qi,  a2,  and  rps  given  by  Theorem  2. 

distributed  to  some  degree  to  all  other  gates.  As  a  con¬ 
sequence,  we  would  like  rp  to  depend  on  both  Ap  and 
some  average  of  all  other  Aps.  We  use  the  geometric 
mean  kAvg  =  ’y/ II  ki  as  an  average  of  the  Aps,  since  it 
has  physical  meaning — it  is  proportional  to  the  theoretical 
minimal  delay  of  the  cycle.  In  this  context,  we  introduce 
the  following 

Theorem  2  For  a  neighborhood  Vp  =  [p  —  p,p  +  rj]  of 
p  >  0,  p  >  0,  and  a  neighborhood  Vk  =  [k  —  r],k  +  p] 
of  k  >  0,  the  values  of  a i,  a2,  Pi  and  P2  that  minimize 
Etn  given  the  of  the  form  defined  by  Equation  13  with 
?p(A;o,  Ap, . . . ,  km-i)  =  j3i  jf-  +  /?2,  where'ii  :  Pi  €  Vp, 
ki  €  Vfc,  and  p  — >  0,  are 

1  1  1 
ai  =  i  m  =  n-j  2  -  ,  and  Si  = 

n  '  m— 1  n  '  m— 1  “ 

Theorem  2  yields  the  optimal  values  of  0:1,02,  Pi  and 
P2  when  the  pi  s  are  in  a  close  neighborhood  of  p,  and  the 
kiS  are  in  a  close  neighborhood  of  k.  We  would  like  to 
verify  now  how  good  these  values  are  in  minimizing  Etn 
when  the  p;s  and  the  kj  s  are  no  longer  close  to  each  other. 
We  use  again  a  numerical  optimizer  to  compute  the  error 
between  the  optimal  and  the  estimated  Etn  for  a  given 
n,  m  and  a  set  of  p,s  and  k^.  We  vary  m  €  [2, 1000], 
n  €  [1, 10]  and  use  three  different  distributions  (uniform, 
uniform-squared,  and  uniform-cubed)  for  p*  €  [1,100] 
and  ki  €  [1,3.3]  (if  we  assume  kni  €  [1,6]  and  kpi  € 
[1,2],  then  with  p,  =  2.5  we  get  ki  €  [6.66,21.95]  or 
equivalently,  using  Property  2,  ki  €  [1,3.3]).  As  for  Equa¬ 
tion  11,  the  observed  errors  are  practically  independent  of 
the  problem  size  to  and  the  distribution  chosen  for  the  p,s 
and  the  fcjs;  the  errors  only  depend  on  n.  Figure  2  shows 
the  relative  error  in  E,  t  and  Etn  for  to  =  31,  n  €  [1, 10] 
andpj  €  [1, 100],  A:*  €  [1,3.3]  chosen  randomly  through 
a  uniform-squared  distribution.  The  average  error  in  E  is 
between  1.7%  and  6.1%,  the  average  error  in  t  is  between 
-3.4%  and  1.7%,  while  the  average  error  in  Etn  is  about 
3.3%  for  n  =  2,  but  increasing  about  linearly  with  n,  ow- 


ingto  the  error  amplifying  artifact  of  Et"  (if  t  =  fo(l+A) 
=4-  t"  ps  tg  (1  +  nA)  for  small  A). 


6.3.  Non-homogeneous  Circuits  (second  form) 

The  main  intended  use  of  Equation  13  in  energy-delay 
efficient  design  is  to  find  approximate  transistor  sizes 
when  n  m  2,  i.e.,  when  voltage  scaling  is  a  design  parame¬ 
ter.  As  Figure  2  shows,  the  equation  stated  by  Theorem  2, 
i.e.,  a  particular  case  of  Equation  13,  does  this  reasonably 
well — i.e.,  within  a  few  percent  of  the  optimum.  On  the 
other  hand,  one  might  want  to  use  such  a  sizing  formula 
for  large  n  as  well — i.e.,  predominantly  delay-only  opti¬ 
mization.  Getting  a  close  approximation  of  Et"  when  n 
is  large  requires  a  very  good  delay  estimate,  since  even  a 
small  error  A  in  t  gets  linearly  amplified  to  nA  in  Et" . 
For  this  reason,  we  study  the  behavior  of  Equation  13  and 
the  delay  estimate  resulting  from  it,  when  n  — >  oo. 

For  now,  consider  a  simpler  problem,  namely  finding 
the  transistor  widths  w ooi  that  minimize  t  given  by  Equa¬ 
tion  7.  This  is  a  special  case  of  the  Etn  optimization  prob¬ 
lem  for  n  — >  oo.  In  [16]  we  have  shown  that  the  optimal 
delay  t0 0  =  mkAvg  is  reached  for  transistor  widths  that 
have  the  property 

Vi  :  =  kAvg  _  (14) 

Woo  i  kj 


We  would  like  the  Wi  s  given  by  Equation  13  to  have  prop 
erty  (14)  for  large  n.  More  precisely, 

tG_|_  1  ^oo(t-t-l) 


lim 

n—>oc 


Wi 


W  o 


(15) 


or  equivalently,  using  ot\  and  ce2  given  by  Theorem  2, 


lim 
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lim  (^0:  ^t  5  •  •  •  i  km— i ) 
n-s-oo  n(k0,  fci,  .  .  .  ,  km-i) 
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(16) 

Condition  16  guarantees  that  the  delay  estimate  re¬ 
sulting  from  Equation  13  is  optimal  for  large  n.  An 
obvious  choice  of  the  r^s  that  fulfills  (16)  is  Vi  : 
ri(ko,  ki,. km- 1)  =  /3wooi ,  where  fi  >  0  is  a  con¬ 
stant  scaling  factor.  The  role  of  8  is  to  normalize  the  Wi  s 
to  the  right  energy  level;  its  optimal  value  is  stated  by  the 
following 


Theorem  3  For  a  neighborhood  Vp  =  [p  —  p,  p  +  r/}  of 
p  >  0,  p  >  0,  and  a  neighborhood  Vk  =  [k  —  p,k  +  rj\ 
of  k  >  0,  the  values  of  a±,  a2,  8  that  minimize  Etn 
given  the  WiS  of  the  form  defined  by  Equation  13  with 
rdk0,  ku....  fcm_i)  =  fiwooi,  where  Vi  :  Pi  €  Vp, 
ki  G  Vfe,  and  r)  — >  0,  are 
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Figure  3.  Accuracy  in  E,  t  and  Et,n  of  the  approximation 
given  by  Theorem  3. 


Assuming  Vi  :  pi  =  p,  and  the  ups  given  by  Theorem  3, 
Equations  6  and  7  yield 

E  =  (1  +  nSifi)E0  and  t  =  (l  + 

where  Eq  is  the  theoretical  minimal  energy  (i.e.  total 
switched  wire  parasitic)  and  too  is  the  theoretical  mini¬ 
mal  delay.  In  [16,  17]  we  have  shown  that  for  a  wide  class 
of  circuits 

Em  (l+n)E0  and  t  m  (l  +  —  ^foo- 

V  n) 

Given  the  value  of  /?  from  Theorem  3,  we  have  that  Vn  > 
0  :  ^7  <  8  <  $2  with  8  —  ^  if  n  ~ >  0  and  8  =  S-2  if 
n  — >  oo.  If  we  choose  8  =  -g -,  the  error  in  E  is  reduced 
by  bringing  E  close  to  (1  +  ti)Eq,  while  if  we  choose 
8  =  S-2,  the  error  in  t  is  reduced  by  bringing  t  close  to 
(1  + 

The  formula  resulting  from  Theorem  3  works  ex¬ 
tremely  well  in  practice  for  small  m ,  i.e.,  it  keeps  the  error 
in  Et"  very  low  for  the  entire  range  of  n,  including  large 
n.  However,  for  m  large  the  accuracy  of  the  formula  de¬ 
teriorates  somewhat  due  to  the  fact  that  E  becomes  con¬ 
sistently  overestimated,  while  the  estimate  in  t  stays  very 
accurate.  This  is  a  consequence  of  the  choice  of  the  rys, 
where  we  have  intentionally  favored  the  accuracy  of  the 
delay  estimation.  For  large  ms,  the  difference  between  -3- 
and  S2  becomes  large  enough  so  that  the  resulting  8  pulls 
E  noticeably  away  from  the  optimum  (1  +  ti)Eq. 

Figure  3  shows  the  relative  error  in  E,  t  and  Et"  for 
the  approximation  given  by  Theorem  3  for  m  =  9  (an  18 
transitions  per  cycle  circuit),  n  €  [1, 10]  andp*  €  [1, 100], 
k,  e  [1, 3.3]  chosen  randomly  through  a  uniform-squared 
distribution.  The  average  error  in  E  is  between  4.4%  and 
6.7%,  the  average  error  in  t  is  between  -0.2%  and  -2.8%, 
and  the  average  error  in  Et"  is  between  1.4%  and  2.3%. 
It  is  interesting  to  point  out  that  for  n  =  100,  the  average 
error  in  E  is  about  1.2%,  the  average  error  in  t  is  about 
-0.003%,  and  the  average  error  in  Et"  is  about  0.5%. 


For  clarity.  Theorems  1,  2,  and  3  were  formulated  to 
refer  to  the  transistor  sizing  problem  of  a  single -cycle  sys¬ 
tem.  However,  these  theorems  can  be  easily  extended  to 
multi-cycle  systems.  We  extend  formula  11,  and  as  a  con¬ 
sequence  Theorem  1,  to  multi -cycle  systems  by  redefin¬ 
ing  PAvg  f°r  each  gate  i  to  be  the  average  parasitic  of 
all  simple  cycles  gate  i  is  part  of.  Theorem  2  extends  to 
multi-cycle  systems  by  substituting  mkAvg  with  too  (the 
minimum  achievable  delay  of  the  circuit).  Given  the  defi¬ 
nition  of  iVroiS  and  PAvg ,  Theorem  3  generalizes  straight¬ 
forwardly  to  multi-cycle  systems,  with  the  only  remark 
that  to — in  the  expression  of  Si  and  S2 — represents  the 
total  number  of  transistors  in  the  considered  circuit,  not 
just  the  ones  on  a  given  cycle. 

Remembering  the  derivation  of  Section  4,  the  values  of 
the  WiS  are  per  gate  i;  but  they  can  be  transformed  into 
the  effective  nFET  and  pFET  sizes  directly,  using  Equa¬ 
tions  3,  4  and  5. 


7.  An  iterative  approach  to  £  fn-optimal 
transistor  sizing 


With  the  help  of  Theorems  2  and  3,  we  can  compute 
approximate  transistor  sizes  of  an  Etn -optimal  circuit.  As 
we  have  seen,  the  approximate  solution  yields  energy  and 
delay  values  within  a  few  percent  of  the  optimum.  How¬ 
ever,  if  the  accuracy  of  such  a  solution  is  not  acceptable 
for  the  given  application,  one  might  wish  to  employ  an  it¬ 
erative  procedure  to  “fine  tune’’  the  initial  transistor  sizes. 

Using  Equation  9,  we  can  compute  Wi — for  a  fixed  i — 
as  a  function  of  the  other  ws.  More  precisely,  if  we  call 
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Figure  4.  Error  in  Et 2  when  exhaustively  simulating  an 
entire  class  of  circuits. 


pre-iteration  value.  This  is  because  =  0,  i.e.,  the 

new  Wi  is  Uf" -optimal  when  all  other  ws  are  fixed  at  their 
current  value.  Secondly,  the  Etn  optimization  problem  is 
convex  in  the  ws,  hence  a  local  minimum  reached  by  the 
iteration  procedure  is  indeed  the  global  minimum. 

To  fully  appreciate  the  benefit  of  the  proposed  itera¬ 
tion  procedure  when  applied  to  the  initial  solution  given 
by  Theorems  2  or  3,  we  exhaustively  analyze  a  particular 
case  of  Equation  8  with  n  =  2,  to  =  5,  Pi  €  {1,2, 3, 4, 5} 
and  ki  €  {1,2,3}.  Figure  4  shows  a  histogram  of  the  rel¬ 
ative  error  in  Et2  between  the  optimal  values  (computed 
with  an  optimization  algorithm)  and  the  estimated  values 
based  on  Theorem  3,  and  also  between  the  optimal  val¬ 
ues  and  the  values  computed  by  one  step  of  the  iteration 
procedure  starting  with  the  approximate  solution  given  by 
Theorem  3.  One  step  of  iteration  assigns  one  new  value  to 
each  Wi.  We  observe  that  the  already  small  maximal  error 
of  the  original  sizing  formula  is  reduced  about  ten-fold  by 
a  single  step  of  the  iteration  procedure.  Of  course,  one  can 
repeat  the  same  procedure  and  get  an  even  smaller  error. 
However,  this  second  step  does  not  have  the  same  impact 
on  reducing  the  error  as  the  first  step  had.  Given  that  the 
transistor  sizes  of  a  real  circuit  are  integer  multiples  of  a 
technology  dependent  constant,  there  is  not  much  point  in 
trying  to  find  the  zero-error  solution.  That  solution  is  un¬ 
likely  to  be  implementable  in  practice,  since  it  will  likely 
have  non-integer  components. 


we  can  compute  Wi  as  the  positive  solution  to  the  cubic 
equation 

+  a2w2  +  a\Wi  +  a0  =  0.  (17) 

(Equation  17  has  a  single  positive  root  for  n  >  1;  this  can 
be  found  using  Cardan’s  method.) 

The  iterative  procedure  starts  with  an  initial  solution 
and  then  repetitively  computes  each  Wi  as  the  positive  so¬ 
lution  of  Equation  17  with  coefficients  computed  from  the 
current  value  of  all  other  ws.  It  is  easy  to  see  that  such 
a  procedure  converges  to  the  Etn  -optimal  solution.  First, 
the  recomputed  value  of  Wi  yields  a  better  Etn  than  the 


We  have  done  several  experiments  in  which  we  tested 
the  dependence  of  the  iteration  procedure  on  the  initial 
starting  point.  We  have  found  that  the  applicability  of  the 
method  strongly  depends  on  the  initial  solution’s  proxim¬ 
ity  to  the  optimal  solution.  Without  a  good  initial  solution 
like  the  one  given  by  Theorem  2  or  3,  the  method  still 
converges  eventually  to  the  optimum.  However,  the  first 
step  of  iteration  yields  a  solution  that  has  an  error  spread 
two  orders  of  magnitude  greater  than  the  solution  result¬ 
ing  from  the  first  step  of  the  iteration  executed  on  the  good 
initial  solution. 


8.  An  algorithm  for  Etn -optimal  sizing 

As  we  have  seen,  the  transistor  sizes  Wi  of  a  system 
optimized  for  Et"  depend  strongly  on  the  wire  parasitics 
Pi.  Unfortunately,  these  parasitics  are  not  known  a  priori , 
since  they  are  attributes  of  wires  that  connect  transistors 
whose  dimensions  have  not  yet  been  found. 

A  two-phase  algorithm  solves  the  problem  of  the  un¬ 
known  parasitics.  In  the  first  phase,  given  the  transistor 
netlist,  each  wire  is  assigned  an  initial  wiring  cost.  The 
more  is  known  about  the  structure  of  the  transistor  netlist 
and  about  a  future  floorplan,  the  more  accurate  such  an 
assignment  will  be.  Based  on  these  initial  wire  parasitics, 
we  can  then  compute  an  initial  estimate  for  the  WiS  with 
the  formulas  established  by  Theorems  2  and  3. 

In  the  second  phase,  we  wire  up  the  pre-sized  transis¬ 
tors  and  extract  the  actual  wire  capacitances  from  the  lay¬ 
out.  With  these  new  parasitics,  we  recompute  the  transis¬ 
tor  widths  Wi.  Finally,  we  may  fine-tune  the  solution  by 
iterating  once  as  described  in  Section  7. 

If  the  accuracy  of  the  final  solution  should  not  be 
deemed  acceptable,  we  can  add  a  pass  through  a  classi¬ 
cal  numerical  optimizer.  Given  the  proximity  of  the  cur¬ 
rent  solution  to  the  optimum,  such  an  optimization  will 
converge  quickly.  In  this  last  phase,  a  more  accurate  tran¬ 
sistor  model  (e.g.,  a  BSIM  model)  can  be  employed,  so  as 
to  bridge  the  gap  between  the  simplified  Uansistor  model 
used  in  this  paper  and  the  actual  Uansistor  behavior. 

9.  Conclusions 

We  have  proposed  a  set  of  analytical  formulas  that 
closely  approximate  the  optimal  Uansistor  sizes  for  cir¬ 
cuits  optimized  for  Etn.  We  have  justified  the  validity  of 
these  formulas  both  mathematically  and  experimentally. 
We  have  proposed  an  iterative  procedure  that  can  further 
improve  the  accuracy  of  the  original  analytical  solution. 
Experiments  show  that,  when  the  procedure  is  applied  on 
the  analytical  solution,  it  converges  much  more  quickly 
then  with  an  arbittary  starting  point.  Based  on  these  re¬ 
sults,  we  have  introduced  a  novel  transistor  sizing  algo¬ 
rithm  for  energy-delay  efficiency. 

10.  Acknowledgments 

The  authors  thank  Catherine  Wong  and  Karl  Papadan- 
tonakis  for  many  stimulating  discussions. 

The  research  reported  in  this  paper  was  sponsored 
by  the  Defense  Advanced  Research  Projects  Agency  and 
monitored  by  the  Air  Force  under  conttact  F29601-00-K- 
0184. 

[1]  Alain  J.  Martin,  Andrew  Lines,  Rajit  Manohar,  Mika 
Nystrom,  Paul  Penzes,  Robert  Southworth  and  Uri  Cum¬ 
mings.  The  Design  of  an  Asynchronous  MIPS  R3000  Micro¬ 
processor.  Proceedings  of  the  17th  Conference  on  Advanced 
Research  in  VLSI,  IEEE  Computer  Society  Press,  p  1 64- 181, 
1997. 


[2]  Alain  J.  Martin,  Towards  an  Energy  Complexity  of  Compu¬ 
tation.  Information  Processing  Letters,  77,  2001. 

[3]  R.  Gonzalez,  M.  Horowitz,  Supply  and  threshold  voltage 
scaling  for  low  power  CMOS.  IEEE  lournal  of  Solid-State 
Circuits,  August  1997. 

[4]  Paul  I.  Penzes,  Alain  J.  Martin,  Energy-Delay  Efficiency  of 
VLSI  Computations.  12th  Great  Lakes  Symposium  on  VLSI, 
New  York,  USA  April  18-19,  2002 

[5]  P.E.Gill,  W.Murray,  M.H.Wright,  Practical  Optimization. 
Academic  Press,  1981. 

[6]  D.P.Marple,  Transistor  size  optimization  in  the  Tailor  layout 
system.  Design  Automation  Conference,  1989,  pp.  43-48 

[7]  I.P.Fishburn,  A.E.Dunlop,  TILOS:  A  posynomial  approach 
to  transistor  sizing.  Proceedings  of  the  1985  International 
Conference  on  Computer-aided  Design,  Nov.  1985,  pp.  326- 
328 

[8]  B.  Hoppe.  G.  Neuendorf,  D.  Schmitt-Landsiedel,  W. 
Specks,  Optimization  of  high-speed  CMOS  logic  circuits 
with  analytical  models  for  signal  delay,  chip  area,  and  dy¬ 
namic  power  dissipation.  IEEE  Transactions  on  Computer- 
Aided  Design,  9(3):236-247,  1990. 

[9]  Y.  Tamiya,  Y.  Matsunaga,  M.  Fujita,  LP  based  Cell  Selection 
with  Constraints  on  Timing,  Area  and  Power  Consumption. 
in  Proc.  ICCAD  Conf.,  Nov.  1994. 

[10]  G.  Chen,  H.  Onodera,  K.  Tamara,  An  Iterative  Gate  Sizing 
Approach  with  Accurate  Delay  Evaluation.  Proc.  IEEE  Int’l. 
Conf.  on  CAD,  1995. 

[11]  Chenming  Hu,  Device  and  Technology  Impact  on  Low 
Power  Electronics.  Low  Power  Design  Methodologies, 
Kluwer  Academic/Plenum  Publishers,  1996 

[12]  M.  Horowitz,  T.  Indermaur,  R.  Gonzalez,  Low -power  digi¬ 
tal  design.  Symposium  on  Low  Power  Electronics,  October 
1994,  pages  8-11. 

[13]  J.  Cong,  C.  K.  Koh,  Simultaneous  Driver  and  Wire  Sizing 
for  Performance  and  Power  Optimization.  IEEE  Trans,  on 
VLSI  Systems,  pp.  408-425,  December  1994. 

[14]  C.  Mead,  L.  Conway,  Introduction  to  VLSI  systems.  Addi¬ 
son  Wesley,  1980 

[15]  J.  Rubenstein,  P.  Penfield.  M.  A.  Horowitz,  Signal  delay 
in  RC  tree  networks.  IEEE  Transactions  on  Computer-aided 
Design  of  Integrated  Circuits  and  Systems  2  (1983),  no.  3, 
202-211. 

[16]  Paul  I.  Penzes,  Energy-delay  Complexity  of  Asynchronous 
Circuits.  Ph.D.  Thesis  (in  preparation),  California  Institute 
of  Technology,  2002. 

[17]  Alain  I.  Martin,  Mika  Nystrom,  Paul  I.  Penzes.  ET2: 
A  Metric  for  Time  and  Energy  Efficiency  of  Computation. 
Power-Aware  Computing,  Kluwer  Academic/Plenum  Pub¬ 
lishers,  2002 

[18]  Paul  1.  Penzes,  Alain  I.  Martin,  Global  and  Local  Prop¬ 
erties  of  Asynchronous  Circuits  Optimized  for  Energy  Ef¬ 
ficiency.  IEEE  Workshop  on  Power  Management  for  Real¬ 
time  and  Embedded  Systems,  Taipei,  Taiwan,  May  29th, 
2001. 


